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Sub-Poissonian light generation in the non-degenerate three-wave mixing is studied numerically 
and analytically within quantum and classical approaches. Husimi Q-functions and their classical 
trajectory simulations are analysed to reveal a special regime corresponding to the time-stable 
sub-Poissonian photocount statistics of the sum-frequency mode. Conditions for observation of this 
regime are discussed. Theoretical predictions of the Fano factor and explanation of the extraordinary 
stabilization of the sub-Poissonian photocount behavior are obtained analytically by applying the 
classical trajectories. Scaling laws for the maximum sub-Poissonian behavior are found. Noise 
suppression levels in the non-degenerate vs degenerate three-wave mixing are discussed on different 
time scales compared to the revival times. It is shown that the non-degenerate conversion offers 
much better stabilization of the suppressed noise in comparison to that of degenerate process. 
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I. INTRODUCTION 

For almost four decades, since the pioneering experi- 
ments of Franken et al. Q and theoretical foundations 
laid down by Bloembergen et al. |0], multiwave mix- 
ing has unceasingly been in the forefront of quantum- 
optical investigations [pM|. In particular, the three- 
wave mixing (TWM) has attracted considerable interest 
as a parametric non-linear process of conversion of two 
sub- frequency (say, wi and UJ2) photons into one sum- 
frequency (uji -\- ijj2 -^ W3) photon, together with the 
inverse process. TWM can be observed in non-linear 
crystals like ADP, KDP, LiNbOs or BaTiOa g. Both 
the total energy, fiuji -\- huj2 = ?ii^3, and momentum, 
?iki -f h'k2 = ?ik3, of interacting photons are conserved. 
TWM is observable for proper orientations of light beam 
polarizations and crystal axes ^ , therefore it can be con- 
sidered as a parametric process. TWM is used for the 
frequency-up conversion if ui —^ 0^3 or the frequency- 
down conversion if ws -^ wi. The process is also use- 
ful for generation of non-classical light such as squeezed, 
sub-Poissonian and antibunched light [p|. 

Before the computer era, the quantum dynamics was 
usually investigated under the short-time approximation 
only. Nowadays, the Taylor series of quantum operators 
can be found for almost any number of terms with the 
help of fast computers and sophisticated software. How- 
ever, these series are usually convergent for short evolu- 
tion times or even for initial time only. Thus, numerical 
quantum methods (see, e.g., ||7|) fail in simulation of the 
long-time quantum evolution. On the other hand, as 
we have shown in |g,^, the method of classical trajecto- 
ries gives very good estimation in the case of strong-field 



interaction (practically, for photon numbers larger than 
10). The computational speed of the classical-trajectory 
method does not depend on the numbers of interacting 
photons and, moreover, for a larger number of photons 
one obtains better precision. Thus, the method is very 
fast and offers a simple substitute for the tedious exact 
quantum numerical calculations. The classical-trajectory 
method enables not only numerical but also some analyt- 
ical predictions, e.g., for stationary Fano factors [||J|] or 
for maximum pump depletion in TWM |10[| . A method 
similar to ours to simulate classical noise in TWM was 
used by Chmela |ll|| . 

In the previous papers, we have studied degenerate 
processes of wave mixing, including the second |8| and 
higher M harmonic generations. Here, we generalize the 
former results for the non-degenerate wave mixing. It 
is well-known that both degenerate and non-degenerate 
TWM can be used for generation of sub-Poissonian light 
||l^,II3]. Nevertheless, theoretical predictions of quan- 
tum parameters, like Fano factor, are most often de- 
rived under the short-interaction (short-time or short- 
length) approximation (see, e.g., O-M), thus valid for 
weak non-linear coupling of the optical fields only. Mo- 
tivated by papers of Nikitin and Masalov [n5| and of 
Bandilla, Drobny and Jex |16 1^, we analyse the long- 
interaction evolution of TWM. The main result of this 
article can be summarized as follows: the TWM can be 
a source of time-stable sub-Poissonian light of the sum- 
frequency mode in the no-energy-transfer regime. The 
deepest noise reduction, with the Fano factor equal to 
5/6, can be observed for the balanced input amplitudes 
ri = r2 = r3/v2. The same degree of photocount noise 
suppression in the sum-frequency mode can be achieved 
for the degenerate TWM. However, the sub-Poissonian 



light produced in non-degenerate TWM is much better 
stabilized compared to that in degenerate TWM. More- 
over, the Fano factors for the sub-frequency modes in the 
non-degenerate TWM are smaller than those for the de- 
generate process. This and other results will be demon- 
strated analytically by applying a method of classical 
trajectories and tested numerically within quantum ap- 
proach. 



O' 



a' 



gt=o 



gt=0.5 



gt=l 



at=1.5 



gt=2.5 



i 






J 




5 




^ 




=; 




5 




• 



-5 


• 



-5 


m 



-5 


m 



-5 


m 



-5 


m 



II. QUANTUM ANALYSIS 

In quantum approach, the non-degenerate three-wave 
mixing can be described by the interaction Hamiltonian 
(e.g., H) 



H = 



hg (aia2al + alalas 



(1) 



where dk and aj, denote, respectively, annihilation and 
creation operators of the sub-frequency (labeled with 
subscripts 1, 2) and sum- frequency (subscript 3) modes; 
g is a non-linear coupling parameter, which is related 
to the quadratic susceptibility tensor x of a given non- 
linear optical crystal and also dependent on the geometry 
of laboratory set-up |j] . 

As in H] , we analyse the quantum Fano factors given by 
-^k = (("I-) ~ {fik)^) / (fik) for a photon-number operator 
fik — d\dk- The light is referred to as sub-Poissonian if 
Fk <1 and super-Poissonian if i^^ > 1 . 

For weak non-linear interactions or short crystal, the 
short-time approximation can be applied for analytical 
predictions of photocount noise suppression with F^ < 1 
p!^ , p^ . The Fano factors are approximated by 

F,,2 = 1 + 2rl {gtf + ^r^r^r^ sinO {gtf + 0{igt)^}, 



i^3 = 1 - ^nr^r^ sine {gtf + 0{{gt)*}, 



(2) 



where rt are the input coherent amplitudes and 9 = 4>i + 
4'2 — i'i is the input phase mismatch. For sm.9 > 0, the 
sub-Poissonian statistics in the sum-frequency mode can 
be observed. For = 0, we find the higher-order short- 
time Fano factor expansion to be 



F:i^l + {rl- Irlrl + Arjrj 



^rlrl 



iaty 



0{{gm. 



It is seen that the sub-Poissonian light in the sum- 
frequency mode is generated for some combinations of in- 
put amplitudes r^ . Since the Fano factors depend weakly 
on time (i.e., in its third or higher-order power), it is dif- 
ficult to observe the sub-Poissonian light generation in 
the short-time regime. 

To analyse the exact quantum dynamics of the TWM 
process beyond the short-time approximation, we apply 
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FIG. 1. Quantum dynamics out of NETR: Cross-sections 
of single-mode Q-functions: (5i(Reai,Imai) and Q2(Rea2, 
Ima2) for sub- frequency modes and (33(Rea3,Ima3) for 
sum-frequency mode at different scaled evolution times for 
initial coherent fields with real amplitudes, afc(O) — rt, set to 
Ti = 6, r2 — 4: and ra = 0. 





gt=o 




gt=0.5 




gt=l 




gt=1.5 




g1=2 




gt=2.5 


3 
-5 


• 


5 

-5 


■^ 




^ 


5 

-5 


f 


5 

-5 


n 


5 

-5 


^ 




-5 5 




-5 5 




-5 5 




-5 5 




-5 5 




-5 5 


5 

a'" 

-5 


• 


5 


-5 


t 




^ 


5 

-5 


r^ 


5 

-5 


e 


5 

-5 


Q 




-5 5 




-5 5 




-5 5 




-5 5 




-5 5 




-5 5 


5 

-5 


• 


5 

-5 


• 




s 


5 

-5 


^ 


5 

-5 


^ 


5 

-5 






FIG. 2. Quantum dynamics in NETR: Cross-sections of 
marginal Q-func tions as in figure 1, but for ri = 6, r2 — 4, 
and rs = rir2/y/rf+r^ ^ 3.328. 
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FIG. 3. Quantum dynamics ottiofNETR: (a) photon num- 
bers (hk) and (b) Fano factors Fk (k = 1,2,3) for coherent 
inputs with real amplitudes n = 6, r2 =4, and rs = 0. 
For longer times, all the modes become stationary with the 
super-Poissonian statistics, Fk > 1. 

the Walls-Barakat method p] of Hamiltonian diagonal- 
ization for the initial coherent states. Quantum analysis 
enables numerical estimation of all statistical properties 
including photocount noise. Complete quantum informa- 
tion of the TWM dynamics can be given by the Husimi 
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FIG. 4. Quantum dynamics m NETR: Fano factors _F\ 
(thin solid), F2 (dashed), and F3 (thick solid curves) ver- 
sus scaled time for different amplitudes of initial coherent 
fields: (a) ri = 6, r2 — 4, and (b) ri — 6, r2 = 2 together 
with rs — Txr^l ^ r\ + r\. Sub-frequency modes become su- 
per-Poissonian, F\^2 > 1- By contrast, the sum-frequency 
mode becomes sub-Poissonian, F-j, < 1. Dotted lines repre- 
sent the classical trajectory predictions of F^' to which all 
the quantum curves tend asymptotically. 



Q-function defined to be 

Q(ai, 02,03) = 7r"''(ai,Q;2,a3|/5|ai,Q;2,a3) 



(3) 



where |ai) (g) \a2) <Xi |a3). In figures 1 and 2, we plot its 
marginal single-mode Husimi Q-functions given by 



Q(Qfe) = / (3(ai, 02,03) Y[ d 



(4) 



where k,m = 1,2,3. The Fano factors, presented in 
figures 3 and 4, were calculated with the help of the 
marginal Q-functions. Due to obvious computational dif- 
ficulties, the exact quantum results can be obtained for 
relatively small numbers (up to few hundreds) of inter- 
acting photons only. 

By analysing the numerical quantum solution we ob- 
serve that the basic features of the photon number evolu- 
tion for the non-degenerate TWM are in agreement with 
those for the harmonic generation processes, as recently 
reported in HJ^. In particular, we observe the so-called 
no-energy-transfer regime (NETR) pS 16 1^, for which 



the energies and intensities of both modes remain con- 
stant in time during the interaction. Although small en- 
ergy flows between the modes appear as a consequence 
of vacuum fluctuations, their influence is negligible for 
strong fields. NETR in the three- wave mixing can be ob- 
served if the amplitudes and phases of the initial coherent 
fields are matched as follows (Eq. (18) in pq]): 



1 

'3 

03 



1 



1 

'2 

P2. 



(5) 



Usually, i.e., for the initial coherent fields not satisfy- 
ing (|5|), all the Fano factors are stabilized in the super- 
Poissonian statistics after a short {gt \rk\ ^ 1) relaxation 
period. Thus, the outputs have high-level photocount 
noise. In figure 1, we present a typical quantum evolution 



of the single-mode Husimi functions Q{ai) and Q{ai) for 
the initial amplitudes ri = 6, r2 = 4 and r^ = 0. The cor- 
responding evolutions of the photon numbers and Fano 
factors are presented in figure 3. 

Different behaviour is observed if the initial phase (/)3 
and amplitude r^ of the sum-frequency mode fulfill the 
condition for NETR. This distinction is clearly seen by 
comparing figures 1 and 2 for Q-functions or figures 3 and 
4 for the Fano factors. In figure 4, the Fano factors are 
calculated for two different pairs of the initial amplitudes 
of sub-frequency modes: (a) ri = 6, r2 = 4 and (b) 
ri = 6, r2 = 2 and the sum- frequency- mode amplitude r^ 
fulfilling (|5|). We observe that all the Fano factor curves 
start at Fk (0) — 1 and after some relaxations become 
stationary at much lower noise levels than those for fields 
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FIG. 5. Non-degenerate vs degenerate TWM oui of NETR: 
Photon numbers (hk) and Fano factors Fk {k — 1,3) are 
obtained for initial amplitudes ri = r2 = 6 and ra = 0, 
whereas (nj.) and F^. are for r'l — 6 and r'^ = 0. All 
quantities with (without) prime correspond to the degener- 
ate (non-degenerate) TWM. 
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FIG. 6. Non-degenerate vs degenerate TWM m NETR: 
(hk) and Fk are calculated for initial amplitudes ri — r2 = 6 
and ra = 6/\/2, whereas (nJ.) and F^. are for ri = 6 and 
r^ = 3. 
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FIG. 7. Revivals and collapses in degenerate TWM vs 
those in non-degenerate TWM out of NETR. Initial condi- 
tions are the same as in figure 5. 
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FIG. 8. Revivals and collapses in degenerate and 
non-degenerate TWM in NETR with the same initial con- 
ditions as those in figure 6. 

out of NETR. The sub-frequency modes remain super- 
Poissonian with F\^2 (t) > 1, whereas the sum-frequency 
mode becomes sub-Poissonian with F3 (i) < 1. The 
most suppressed noise is observed for the balanced in- 
puts, given by ri = r-i and rs — r\l\/2, as a special case 
of condition (H). For those inputs, the Fano factor in the 
time limit tends to F3 (t -^ 00) w 0.83. 

As we have shown in Refs. |8|j^, the same degree of 
the Fano factor can be obtained in degenerate TWM. 
Thus, one can address the following objection: What is 
the point to study the same parameter in the closely re- 
lated non-degenerate version? First, we stress that the 
same Fano factor is obtained in a special case only: for 
the sum-frequency mode in NETR for long-interaction 
times and high-intensity fields. By contrast, these fac- 
tors are distinct for the sub-frequency modes under the 



same initial conditions. In general, the results even for 
the sum-frequency mode in these two processes are dif- 
ferent for initial conditions either not fulfilling (||) or for 
lower intensity fields or different time periods. Second, 
we will show that the Fano factors for non-degenerate 
TWM are better stabilized than those for degenerate 
TWM for much longer evolution times. This is an im- 
portant advantage of the non-degenerate conversion. 

For better comparison, let us analyse in detail the de- 
generate TWM described by 



H' ^ 



hg' [alal 



1 -t2- 



(6) 



For clarity, observables calculated for the degenerate 
TWM are marked with prime to distinguish them from 
those for model (ll]) and we keep subscript '3' (not '2') for 
the sum-frequency mode. Hamiltonian (|6|) formally dif- 
fers from (m) in the assumption of ai = 0,2 only. But we 
also put g' — (?/v2 for better synchronization of oscilla- 
tions in photon numbers (nj^,) and {fik)- Sub-Poissonian 
statistics in degenerate TWM was analysed in our for- 
mer works MM- In figures 5 and 6, we compare exact 
quantum evolutions of the mean photon numbers and 
Fano factors for degenerate {{n'l^) and F^. for k = 1,3) 
and non-degenerate {{fik) and F^) TWM for the same 
initial conditions. We observe similar behavior for short 
times gt <C 1 only. For longer times (also for gt < 1) 
results arc not equal by no means. Different predictions 
of quantum evolutions for models (0) and (O) come from 
different commutation relations: [01,02] = for non- 
degenerate TWM and [01,03] = [01,0]^] = 1 for degener- 
ate TWM. As a result, constants of motion are different: 
^totai(*) = "U*) + 2n(j(i)=const for degenerate TWM, 
whereas hi{t) — n2{t)=const and ni{t) + n3(i)=const 
for non-degenerate TWM. We note that Ntota\{t) = 
hi{t) + fi2{t) + fisit) ^const in the latter process. Evo- 
lutions, presented in figures 5 and 6 for degenerate and 
non-degenerate TWM, are distinct in amplitudes and fre- 
quency of oscillations as well as in the level of their "sat- 
uration" . For example, the limiting value of the sum- 
frequency Fano factor for degenerate TWM is lower than 
that for non-degenerate TWM for evolution out of NETR 
(see figure 5(d)). While the sub- frequency Fano factor in 
NETR is considerably higher for degenerate compared to 
non-degenerate TWM (see figure 6(c)). 

In figures 1-6, we have analysed the time regime which 
is long compared to the typical interaction times for 
known crystal lengths. However, the time is short com- 
pared to the revival times for such systems. The ques- 
tion arises about the photon-number noise suppression 
on such a long scale. This analysis will show an advan- 
tage of non-degenerate over degenerate TWM related to 
the stabilization of the suppressed photon-number noise. 
In figures 7 and 8, we present the long-time evolution 
for < gt < 100 of the mean photon numbers and Fano 
factors for both non-degenerate {{fik) and Fk) and degen- 
erate ((n'j.) and Fl.) TWM. We observe that the revivals 



are strongly pronounced for (i) degenerate TWM (right 
figures 7 and 8) compared to non-degenerate process (left 
figures), and (ii) outside NETR (figure 7) rather than 
in NETR (figure 8). Thus, the non-degenerate TWM in 
NETR exhibits the highest stability. Even for longer evo- 
lution times as 100 < gt < 1000 of non-degenerate TWM, 
the oscillations are similar to those for < gt < 100 and 
it is hardly difficult to classify them as a typical revival. 
Oscillations in (hk) are of order 10"'^ and in Fk are of 
order 10~^ even for such small initial amplitudes equal 
to «! = a2 = 6 and a^ = 6/-\/2- Our analysis is re- 
stricted to initial coherent inputs. It is worth noting 
that the revivals are much stronger for quantum input 
fields like, e.g., Fock states. In conclusion, although the 
degenerate and non-degenerate TWM lead to approxi- 
mately the same photon-number noise suppression in the 
sum-frequency mode for NETR (see figure 6(d)), the non- 
degenerate process offers much better stabilization of the 
suppressed noise for long evolution times (compare fig- 
ures 8(c) and 8(d)). 

In the next sections, we will apply an approximate 
method of classical trajectories to explain the extraordi- 
nary stabilization of the observed photocount noise and 
to estimate analytically the level of noise suppression for 
NETR. 



III. CLASSICAL ANALYSIS 

Complete quantum solution of the model given by 
Hamiltonian (w) can be found numerically only. Yet, in 
a special case for strong fields, analytical results can be 
obtained by applying approximate classical methods. 

In analogy to quantum Hamiltonian (Q), the classical 
model of non-degenerate TWM can be described by B : 



H = g {aia2Ctl + c.c. 



(7) 



where ak are the complex amplitudes of the fc-th mode 
and g is a non- linear coupling parameter. From (R), one 
readily obtains the following complex differential equa- 
tions 



ai = -iga2a3, 
^2 = -igalas, 
<^3 — ~igaia2- 



(8) 



It is easy to show by comparing (||) with equations (12) 
and (13) from Ref. Q that the classical models for the 
degenerate and non-degenerate TWM are equivalent for 
ai = a2 and arbitrary evolution times. To get equations 
of motion for degenerate TWM, it is enough to replace 
g by V^g' and a^ by V^a'^ in (||). As was discussed in 
the former section, the quantum evolutions of degenerate 
and non-degenerate TWM are equivalent for gt <C 1 only. 
By introducing real amplitudes and phases, a^ = 
TkC^'^'' , equations (0) can be transformed into the follow- 
ing four real equations 



ri 
r2 



■ -gr2r3sm0, 
-grirssine, 
grir2 sin0, 
Vir2 rirs 



9 



rs 



r2 



r2r3 
ri 



cost 



(9) 



where 61 = (/)i -I- 02 — 03 is the phase mismatch. The 
system (P) has three integrals of motion 



El 

E2 

K 



r^ + Tg = ni 

2 I 2 
r2 + Tg = 722 

rir2rj, cos 9. 



-n3, 

(10) 
equation for the re- 



By extracting ri , r2 and 9 from 
maining amplitude r^ reads as 

{rsh/gf + K^^ rl {El - rj) {E2 - rj) (11) 

or, equivalently, for the intensity n^ — r| as 

ih3/2gf = n3 {El - ng) (^2 - "3) - K^ 

= {a - ns) {b - ns) {n3 - c) , (12) 

where the numbers a > b > c are the roots of cubic 
equation 713 {Ei — n^) (£'2 — n^) — K^ = satisfying the 
conditions abc = K^ , a + b -\- c = Ei -\- E2, and ab + ac + 
be — EiE2- Then, the solution for n^ (t) is found to be 



n^ {t) = c+ {b — c) sn [Va — cgt + 0o, k\ , (13) 

where sn(w, fc) is the Jacobi elliptic function with k = 
\/ ^^ and 00 is the initial phase given by the elliptic 
integral of the first kind 

dx 



(bo^F{z,k) 



Vl-fc2 



• 2 

sm X 



(14) 



where z = arcsin 1/(713(0) — c)/(6 — c). One observes 
that 71,3 is a periodic function oscillating between the val- 
ues c and b with the period given by T = 4F {^,k) /g. 

In two special cases, solution (13) reduces to the ele- 
mentary solutions: 



713 (t) = f tanh (rgt) 
for ri = r2 = r and r^ = 0, and 

"-3 (t) = r'^sech^{rgt) 



(15) 



(16) 



for ri = r2 = and ra = r. Another elementary solution 
is obtained for the initial fields fulfilling conditions (^). 
In this case, the solution reads as {k = 1,2, 3) 



ak (t) = Tk exp 



.rir2r3 



(17) 



which describes the classical no- energy-transfer regime 
(NETR) Q , since the amplitudes and energies of all in- 
teracting modes remain constant, i.e., nk(t) = \ak{t)\'^ = 
rl- We conclude that NETR observed in the quan- 
tum numerical analysis presented in former section cor- 
responds to the classical solution (|l7|). 



IV. CLASSICAL TRAJECTORY ANALYSIS 

Classical solutions, as presented in the preceding sec- 
tion, do not describe quantum noise. Nevertheless, they 
can be used for simulation of quantum noise if the ini- 
tial complex amplitudes are chosen randomly. This ap- 
proach, referred to as the method of classical trajectories, 
has been applied successfully in a description of noise 
in various quantum-optical phenomena p- [ll| , [l5| [20|] . By 
analysing Q-functions and Fano factors, we will show 
that the method of classical trajectories properly sim- 
ulates photocount noise in the TWM processes. 

To calculate statistical moments, like the Fano factors, 
one needs to analyse the classical evolution of each pro- 
cess (trajectory) separately and then to average the mo- 
ments over all the obtained trajectories. The classical 
Fano factor, defined to be 



F — 



(18) 



can be obtained by the classical trajectory averaging. We 
denote this averaging by bar to distinguish it from quan- 
tum ensemble averaging indicated by brackets. We will 
apply the method of classical trajectories along the lines 
of the analysis presented in Q . We choose the initial am- 
plitudes to be ak = Tk and blur them with the Gaussian 
noise, which results in 



ak = rk +Xk +iy/c, 



(19) 



where X)~ and y^ are mutually independent real Gaussian 
stochastic quantities with the identical variances a^ = 
1/4. We assume further that the fields are strong, i.e., 

Tk > 1. 

By substituting (|lj) into (^0|) , the integrals of motion 
can be expressed as 



El = |ai 



laal^ =rl+rl + ci, 



^2 = |a2|' + |a3|'=r-2-Kr2+c2, 
K — Ke {aia2a*^) = ?'ir2r3 + di, 



(20) 



where ci, C2 and di are small corrections of E'l, E2 and K , 
respectively, depending on r^, Xk and yk- To eliminate 
the linear term in the RHS of (E2[), we substitute 



1 



ny. 



1 



1 



-E2 + -El - -\ {Ei - E1E2 + Ei) + e 



n?,a 



(21) 



where e is a small correction of stationary value 7130 + b. 
Under the strong- field approximation (rk 3> 1), one can 
neglect the small correction e'^ and the RHS of ( p^ ) can 
be approximated by the quadratic function fi^ (a^"^ e^)- 
Thus, one gets simple equation 

e' = {2gnf {a^ - e^) , (22) 

which leads to the solution of (|lj) in the following form 
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FIG. 9. Classical simulation of typical quantum dynamics 
out of NETR; Clouds of 10,000 points representing marginal 
Q-functions for the same initial conditions and times as in 
figure 1. 
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FIG. 10. Classical simulation of quantum dynamics in 
NETR for the same cases as in figure 2. 



na (t) = ?i3o + b + a sin (2gQt + tp) 



(23) 



where n^o 



rir2/{ri -l- r|). The coefficients a, b 



and fi, together with ci and C2, are complex functions 
of Tfe and noise parameters Xk and yk- With the help 
of integrals of motions, given by (10|), solutions for other 
modes (fc = 1,2) can readily be found as 



rik (t) — Ek — n3{t) — rl + Ck — b — a sm{2gQt + ip) . 



(24) 



We observe that all the three solutions, given by (|23D and 
([2^), are of the form of large constants slightly perturbed 
by the same harmonic function. In figures 9 and 10, we 
present classical simulation of quantum dynamics by cal- 
culating time evolutions of 10,000 points in phase space 
according to classical equations of motion. These rep- 
resentations correspond to the Husimi Q-functions pre- 
sented in figures 1 and 2, respectively. By comparing 
figures 1 and 2 or, equivalently, figures 9 and 10, we ob- 
serve two distinct types of evolution determined by the 
initial amplitudes to be in NETR or out of it. 

The classical and quantum descriptions are principally 
different in detail. Thus, our plots of the Q-function 



based on the exact quantum solution of TWM (figures 
1 and 2) and those obtained by an approximate classical 
simulation (figures 9 and 10) also differ in detail. The 
discrepancies are more pronounced for lower amplitude 
inputs and longer interactions. Moreover, the methods 
of graphical representations are different: a topograph- 
ical picture of Qfe(a) versus a cloud of classical points. 
Nevertheless, "it is surprising how close the clouds of dots 
are to the Q-function" |13] . The clear correspondence be- 
tween figures 2 and 10 or 1 and 9 justifies our application 
of the classical trajectory approximation. 

As the next step of the classical trajectory method, 
one has to perform averaging of solutions (^ and ( p4[ ) 
to calculate the desired moments. We find that the mean 
values of the parameters occurring in solution (|24| ) are 
6 = cT = C2 = and 



n^^lrl+rl 



/r| + r'lr2 + r\ 



(25) 



together with their mean quadratic moments 



P^^ 



2A 



(ci - bf = ^-(4r° + lQr\ri, + llr^r^ + 7rfr° + 2r=) 



{c2-hr 



'2 

2A 



4„4 



„2„6 



{2rl + 7r\r^2 + 11^1^2 + 10^1^2 + 4^: 



2h 



given in terms of the auxiliary function 



^2^2 



A 



8(rt + rfr2^+r|) 



4\2 



(26) 



(27) 



The phase (/s can be obtained from ( |23[) at i = 0. Thus, 
the photon-number mean values are simply equal to fifc — 
r\ and their variances are given by (fc = 1, 2, 3) 



{cu-hf^\a^ 



in terms of the moments (Eq) and C3 = 0. The 
term sin^ {2g^t + ip) has simply been estimated as ^ 
for sufficiently large t, when fik and Fk become time- 
independent. Thus, we arrive at the following Fano fac- 
tors 






cl _ 1 o /l/'„2 I „2\2 



F^' ^ 1 - SAiri + riY 



(29) 



where A is given by (^7|). As one could expect, the for- 
mulas for F^:^ are symmetric with respect to exchange 
of the subscripts 1 < — > 2. We finally conclude that the 
TWM in the no-energy-transfer regime can be a source 
of time-stable sub-Poissonian light in the sum-frequency 
mode as described by the Fano factor 




4 6 8 10 

ri/r2 

FIG. 11. Classical predictions of Fano factors F^' 
(k = 1, 2, 3) versus ratio ri/r2 of the input coherent-field am- 
plitudes. Small circles represent the quantum Fano factors 
Fk obtained from the exact quantum solutions presented in 
figures 4 and 6. 



^3 



(P) 



3P(1 + P)^ 

(l+p + /92)2 



< 1, 



(30) 



depending on ratio of the mean intensities of initial co- 
herent fields defined by p = '^i/'"! if J'l > or p = r2/rl 
if r2 > 0. The sub- frequency fields become super- 
Poissonian with the Fano factors 



Ffip) = F^\l/p) = 1 + 



p(5 + 5p + 8/j^) 
8(l + p + p2)2 



> 1. 



(31) 



In figure 11, the classical predictions of the Fano factors 
are depicted as a function rxjri- By analysing (pfl) and 
figures 4, 6, and 11, we conclude that the sum-frequency 
mode solely is sub-Poissonian (F3' < 1) and the strongest 
noise suppression is obtained for r\= r2= r^/y/2, when 
-^3' — 5/6 ~ 0. 833. For highly unbalanced input in- 
tensities ri <C ^2 or ri ^ r2, all the Fano factors ap- 
proach unity F{^' « Fj' ~ Fg' — > 1. Mutually equal 
Fano factors, estimated by Ff = Ff = 5/4 = 1.25, 
are predicted for the balanced inputs ri — r2 = r^^j \f2. 
The maximum values of the Fano factors, estimated by 



(28) Ffi 



Ff 



max = 1. 255, are obtained for slightly un- 



balanced inputs: r\ — 1. 136r2 to maximize F^^ and for 
ri = 0. 881r2 to maximize F|'. 

We have predicted in [g|,|| the stationary sub- 
Poissonian Fano factors for the second (F^' — 5/6) and 
third (F|' — 13/16) harmonic generations within NETR. 
The minimum value of the sum- frequency- mode Fano fac- 
tor for the non-degenerate TWM is the same as that ob- 
tained for degenerate TWM (i.e., second-harmonic gener- 
ation) 0] , but higher than that for degenerate four- wave 
mixing (i.e., third-harmonic generation). However, for 
the sub-frequency modes, the Fano factors for the non- 
degenerate TWM are smaller than those for the degener- 
ate cases, namely F{^2(1) = 5/4 instead of 3/2 and 29/16, 
respectively. 

In figures 3 and 4, we have compared evolutions of 
the exact quantum Fano factors Fk (depicted by solid 
or dashed curves) with their classical estimations, F^' 



(dotted lines). One observes that all the curves start 
at i^fc = 1 for initial coherent fields and after some re- 
laxations become quasi-stationary with Fk ~ F^', given 
by (|2^). It is worth noting that very good estimation 
is achieved even for relatively small amplitudes, e.g., 
Tk < Q- We conclude that the conditions for NETR 
in quantum dynamics and suppression of the observed 
quantum noise levels are well explained by the classical 
trajectory method. 

Finally, we will compare scaling properties of the Fano 
factors in their dependence on light intensity and initial 
amplitude for the degenerate and non-degenerate TWM. 
Drobny et al. |l^ calculated the scaling laws under the 
truncated Wigner approximation for the maximum sub- 
Poissonian photon-number noise in TWM. Their for- 
mulas are valid also in the limit of r = ai(0) -^ oo. 
Here, we focus on nonlinear fits for finite ranges of r 
only. Let us investigate the scaling properties of the 
maximum sub-Poissonian character of the sum-frequency 
mode corresponding to the first minimum of the F^ and 
Fg curves in figures 3(b), 4, 5(d), and 6(d). In fig- 
ure 12, we plot the exact quantum numerical values of 
m.intF3{x,t) = F^ix^taun) or mint F^{x,t) = F^(x,t'^-^J 
as a function of initial amplitudes x = ri = r2 = r[ and 
of intensities x = (n3(imin)) = ("■3(^min))- We fit those 
minima with the exponent and polynomial functions of 
parameters listed in tables I and II. The errors, given in 
the last column, are estimated by the standard deviation. 
We observe that the scaling laws ax'' for non-degenerate 

TABLE I. Scaling laws ax'' for non-degenerate TWM 
(mintF3(2;)) and degenerate TWM (mintF3(a;)) out of 
NETR. Initial conditions are 1 < ri = r2 < 10, 1 < rj < 20, 
and r-3 = 7-3 = 0. 



No. 


fitted function 


a 


b 


error 


1 


mintF3(ri) 


0.8819 


-0.1254 


0.0004 


2 


mint Fa ((ns)) 


0.8560 


-0.0572 


0.0001 


3 


mintF3(ri) 


0.7694 


-0.0906 


0.0003 


4 


mintF:^((n^)) 


0.7352 


-0.0427 


0.0003 
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TABLE II. Polynomial fit {ax'^ + bx + c)/x^ for 
non-degenerate (mint F3(2;)) and degenerate (mintF3(a;)) 
TWM m NETR. Initial conditions are 1 < ri = r2 < 10, 
l<r[< 20, rs = rl/^/2, and r'g = r[/2. 
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FIG. 12. Maximum photon-number noise suppression for 
non-degenerate (curves A and C) and degenerate (B and D) 
TWM: Time-minimized Fano factors, mint F3 and mint F3, as 
a function of (a) initial amplitudes r = ri = r'l and (b) mean 
photon numbers (hs) — (n'j). Curves: A (with diamonds) 
for initial amplitude rs = 0, B (with circles) for r'j = 0, 
C (with solid triangles) for rs = ri/^/2 — r/\/2, D (with 
empty triangles) for r'3 = ri/2 = r/2. The marked points 
are obtained from the exact quantum solutions and are fitted 
with the functions given in tables I and II. 

and degenerate TWM give good approximation of the 
exact values out of NETR only. However, the ax'' law 
fails to describe with good precision mint F3 (x) and 
mintF3(a;) in NETR at least for intensities up to 100 
photons. Thus instead of the exponent law we apply the 
(inverse) polynomial fit of the form {ax^ + bx + c)/x^, 
where we introduce x~^ in relation to the definition of 
Fano factor. By contrast to the exponent fits, the poly- 
nomial laws give very good predictions of the maximum 
sub-Poissonian behavior at least for 1 < ri = r2 < 10 
and 1 < rj^ < 20 as seen in figure 12 and table II. The 
scaling laws of Drobny et al. [|^ differ slightly from ours 
presented in table I. The minor differences in the fitted 
parameters result from different ranges of r used in the 
fitting procedures and from application of the truncated 
Wigner approximation in Ref. |17| ] compared to our exact 
quantum method. 

On the other hand, the Fano factors for the balanced 
(ri = r2) non-degenerate TWM under NETR conditions 
for long times and high intensity fields do not depend 
on light intensity, which follows from equations (30) and 
(pn|). Similarly, there are no scaling properties of the 
Fano factors in the degenerate TWM for long times and 
high intensity fields in NETR, i.e., under the same condi- 
tions as those assumed in our classical trajectory analysis. 



No. 


fitted function 


ri 


a 


b 


c 


error 


1 


mint F3(ri) 


> 1 


0.5474 


0.1104 


0.2956 


0.0003 






> 1 


0.5519 


0.0643 


0.3894 


0.0007 


2 


mintF3({n3)) 


> 1 


0.5562 


0.3143 


-0.1618 


0.0001 






> 1 


0.5559 


0.3240 


-0.2057 


0.0008 


3 


mintF3(ri) 


> 1 


0.5495 


0.1195 


0.3510 


0.0003 






> 1 


0.5541 


0.0406 


0.5739 


0.0012 


4 


mintF^({n^)) 


> 1 


0.5558 


0.2057 


-0.0607 


0.0001 






> 1 


0.5556 


0.2121 


-0.0783 


0.0005 



V. CONCLUSIONS 

We have analysed the long-time interactions in the 
non-degenerate three-wave mixing. To the best of our 
knowledge, our quantum analysis is the first presenta- 
tion of the exact and completely quantum solution of 
the nondegenerate three-wave mixing. In literature, a 
special solution can be found for initial sub-frequency 
fields with zero amplitudes ai ~ a2 = only. The 
no-energy-transfer regime for proper choices of ampli- 
tudes and phases of the initial coherent fields has been 



observed. We have compared the evolutions of the 
Husimi Q-functions and their classical trajectory simula- 
tions for processes in the no-energy-transfer regime and 
out of it. We have shown numerically in the quantum- 
mechanical approach that the three-wave mixing in the 
no-energy-transfer regime exhibits the time-stable pho- 
tocount statistics. This phenomenon was explained ana- 
lytically by applying the method of classical trajectories. 
We have shown that the sub-frequency modes become 
super-Poissonian with the Fano factor Fi_2 > 1, whereas 
the sum-frequency mode becomes sub-Poissonian with 
i^3 < 1. We have found that the most suppressed 
photocount noise, given by F^ « 5/6, is obtained for 



the balanced initial intensities 



of the sub- 



frequency modes and the sum-frequency intensity equal 
to r| = r1/2 as determined from condition d|) for the 
no-energy-transfer regime. Scaling laws and polynomial 
fits for the maximum sub-Poissonian behavior have been 
found for different processes and initial conditions. We 
have compared in detail the non-degenerate and degener- 
ate conversions on time scales short and long compared 
to the revival times. We have observed that the non- 
degenerate three-wave mixing, contrary to the degener- 
ate conversion, exhibits stabilization of the suppressed 
photon-number noise even on the revival time scale. 
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